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Abstract 

In this paper, we demonstrate the existence and uniqueness a semi- 
analytical solution of an inverse heat conduction problem (IHCP) in the form 
> Up = Uge in the domain D = {(2,t)|}0<a<1,0<t< TH, u(z,T) = 
f(x), u(0, t) = g(t), and uz(0,t) = p(t), for any 0 << t < T. Some numerical 
experiments are given in the final section. 
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1 Introduction 


The procedure to solve an IHCP is very important in determining unknown tem- 


perature histories and heat flux from known values in the body, which are usually 
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measured as a function of space and time. Especially, a direct measurement of 
the heat flux or temperature in a boundary or initial time of hot body is almost 
impossible. Therefore, recent studies of IHCPs have been numerically treated and 
extended of multiple dimensions with the help of computing architecture. Some 
numerical and theoretical approaches to IHCPs are summarized in [4], [2]. Such 
a procedure using an exact solution are given by Burggraf [3]. It has been shown 
that, if an error is made in known boundary condition, then there will be some 
errors in unknown heat flux of other boundary. A lower bound of this error can 

1 


be estimated by FAT sinh( > 


vation that small values of time At can produce large error in surface flux. In 


). These results are consistent with earlier obser- 


this paper, we apply a finite difference method of semi-implicit type for ou and 
use a parameter @, for driving a stable and convergent solution to the IHCPs. 


Now, suppose that for any given t,0 <t< T, u(z,t) € C*[0, 1] and satisfying 


ACT) Ste lat) nD =O ea SD Or, (1.1) 
Ms = fa). OSes 1, (1.2) 
u(0,t) = g(t), Oe eer ie (13) 
(OE) ple): <= TT (1.4) 
biel Vet) ee TCG) Os £07, (1.5) 
ula): -= oO), O<a2 <1, (1.6) 


where f(x), g(t), and p(t) are piecewise-continuous known functions, T is a 
given positive constant number and h(t), (x), and u(x,t) are unknown func- 
tions, which remain to be determined. 

In the next section, we discrete the variable t and reduce (1.1) — (1.4) toa 
system of linear, nonhomogenous second order differential equations. Stability 
and convergency of this method is studied in section 3. Some numerical result 


and discussion are given in section 4. 
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2 A Numerical Solution 


In this section, we discrete the variable t to approximate the solution of (1.1) — 


(1.4). 

Let MEN, Aty = aa and t; = iAt, fori = 0, 1, ..., M. For the solution 
u, we define u;(z) = u(z,t;), fori =0, 1, ..., MW. Similarly for a given sequence 
of functions {u;(x)| i= 0, 1, ..., MZ} we use t;(x) instead of the approximate of 


u;(x). Putting 


tigi (@) = 0;(x) + (Om OHNE; ts) 


ot 


O(a, ti41) 
Ot 


(Ou — 1) )Atu, (2.1) 


for 014 > 0, then by using (2.1) into (1.1) — (1.4) we obtain a system of linear 
nonhomogenous second order differential equations with given initial conditions 


in the form 


Ou Atty (x) + G(r) = tii (x) + (Ou — 1)Atyiy, 1 (2), (2.2) 
ti(0) = g(ti) = gi, (2.3) 
(0) = p(ti) = pi, (2.4) 


fori=0, 1, ..., M—1, where to(x) (0 < x < 1) and &,(1),i=0, 1, .., M-1 
are unknown. 
Clearly ti (x) = uys(x) = f(x). Using these assumptions, the problem (2.2)— 


(2.4) has a solution of the following form 


x a(1) — V/(1) x 
a(x) = g; cos 4+ + ; sin + Fi(x), (2.5 
ie VouAtu W!(1) AS ay AN 1222) 
fori =0, 1, .... WM —1, where 
1 a Lr-—s 
F,; = —_____ th; in —————— d 
©) = a faaaig J, Hulse a 
=}, es Atm pi41 sin Z om — 1s (x) 
ie M Pi+1 Guts 6 i+1 
Oy —1 x , 
ES SS A Aes a 2.6 
Ou Ji+1 COS i Bia z (2.6) 


V; and W; are the solutions of the following problems, respectively, 


Oy AtyV;" (x) +Vi(z) = ti41(x) + (Ou — DAty tt, (2), 
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and 


WO) =] 1. 40). town =. (2.8) 
Clearly, the heat flux t,(1,¢;) may be obtained from (2.6) of the form 
di2(1, ti) = pi Wi (1) + Vi (1), i=0,1,..., M-1. (2.9) 


Now, for each n € N, if @yAty 4 (nz)~? and f” (zx) is a piecewise-continuous 
function in [0,1], then the system of solutions (2.5) are unique [1]. 


The above result may be summarized in the following statement. 


Theorem 2.1. If, for each n EN, OyAty # (nt)~?, and f" is a piecewise- 
continuous function in [0,1], then the system of differential equations (2.2)-(2.4) 
has a unique solution. 


Proof. See the analysis preceding the of above theorem. 


3 Stability and convergency of solution 


In the next theorem, the convergency of the solution (2.5) to the unique solution 


will be shown. 


Theorem 3.1. Jf O47, Atyand f satisfy the assumptions of Theorem 2.1 and 


2 
|e ee | <C < oo for any t € [0,1], where C is a positive constant number and 


Ou = Bm At ™ such that By > 0, 1+ nae <ay <1 for any M > 3, then 


the solution of the system (2.2)-(2.4) is convergent to the unique solution of the 
problem (1.1)-(1.4), for anyO<a< 1. 
Proof. For alli=0, 1, .... M—1, we have 


OuAty ujy(%) + u(x) =  ugzi() + (Ou — 1)Atyy uj, 1 (2) 
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1 Pu(x, &) u(x, mi) 
= 5 Ati (On are (Ou 1)? ate ), (3.1) 


where t; < & < t+ Oy Atyy and tiga << tiga + (Ou — 1)Aty. 


Now, if we put 
6(e) =e). tle) for t= 0; Tp og My, (3.2) 


then e;(x) satisfies. 


Oy Atu et! (x) + e;(x) €i41(2) + (Ou _- 1)Atyy e¥41(z) 


1 2 2 Pula, &) 2 u(x, m) 
_ pAtu (Ou Ore (Ou 1) Or ); (3.3) 
and e;(0) = e,(0) = 0, from which, we conclude that 
1 cs . L-— 8s Ou —1 
e; (2) mea €j41(s) sin Aah 8 Oxi €441(2) 

2. tyes ae 2 0'u(s, &) 20°u(s, mi), _. LS 

5 om ine F (is aR (Ou — 1) PR ) sin aah 
= Fea) + I;2() + 1 (2) 0, Ls xees M-1. (3.4) 


Clearly, the integrand in [;,3(a) denotes the truncation error and the other 
terms J;1(a) and I;,2(x) show that, errors of initial and boundary data for the 
problem (1.1) how to propagate. In remaining of the proof, we consider two cases: 


Case I. If 


e;(0)=e)(0) = 0, i=0,1,.., M—-1, 
eye = 0, Ode ae ap (3.5) 
then clearly one may conclude that 
[Li,3(2)| < C(@uAty)?”, for i=0, 1, .., M—1, (3.6) 


where C is defined in Theorem 2.2. Thus, ey—1(%) = In—1,3(a) and one may 


show that 


1 
leara(o)| S$ (S057 At +1 - OA lenr—i-a(@)| + asia) 
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i-1 


- ior a2 
= we —6y + 5 Ont Atay) r—i-e.a(@)| 
k=0 


i-1 
1 
< C(OuAtu) SS Nsoy Aty +1—- Oy)", 4=1, 2. M. (3.7) 
k=0 


Clearly, for any fixed value 0y,, limit |e;s_;(x)| will not be zero, when Atyy 
is vanished, which implies the divergency of the solution. If we choose Oy, = 
Bu Ati, ™ for 0 < ay <1 and By > 0, then for each M > 3, and ay which 


satisfies in the following inequality 


In(2G,,r) 
1 4+ ———_ < < 
= In Aty Som $1, 
we obtain that 
lenr—s(x)| < PAB RO em) (3.8) 
where 
i-l 
=CS (1-65; + SOP Mts} Dyes 
k=0 


Consequently, if Atay tends to zero, then |e,s_;(x)| is vanished and w(x) for 
i=1, 2, ..., M, convergences to the exact unique solution of the problem (1.1)- 
(1.4). 


Case 2. In this case, let us suppose 


G = Gt+E&,;, 
Di = Pi + G2, for 1=0, 1, teey M 1, 


and f(x) = f(x) + €(x), then by using (2.5), (2.6), (2.7) and (2.8), we obtain 


E(x) = € 1 CoS Ts + 6,27 @uAty sin Ts 
1 % r-s 
a si Atm €+1,2 Sin Vou 
_ Gar V Atm €+41,1 COS Veuaun + I,3(x), (3.9) 
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where e;(z) is a global error for i= 0, 1, .... M, and ey(x) = €(2). 
Now if |e(x)|, Je1], and |é;2| tend to zero for each i= 0, 1, ..., M, then 
i-1 1 
lenr—s(a)| < $0 (1 — On + 5 OM Ati)” [Im—i+k,3(2)|, 
k=0 
and 


[I—itk,3(2)| < Ce Asian) AiG: 


Finally, |e,_;(xz)| vanishes for alli =0, 1, ..., M. @ 


4 Numerical Examples 


In this section we will present simulated cases to evaluate the capability of the 


proposed robust input estimation scheme. 


Example 4.1. Assume that 


FG) Sa, 
g(t) = 2t, 
p(t) = 2, 

[ ae 


Obviously, u(a,t) = 2? + 2t is an exact solution of the problem. Now, we use 
our numerical method to this problem. For « = 1, Atjyy = 0.02, agg = 0.8, and 


Bur = 0.9, the result are given in the following table. 


t | @Q,t) | u(,¢) Sainte lie ug (1,t) | &e(1, 4) el 
0 | 0.993614] 1 0.0063 2 | 1.98329 0.0083 

0.2 | 1.38377 | 1.4 0.019 2 | 1.96382 0.018 

0.4] 1.77821 | 1.8 0.012 2 | 1.95272 0.023 

0.6 | 2.1803 | 22 0.0089 2 | 1.95764 0.021 

0.8 | 2.58875 | 2.6 0.0043 2 | 1.97603 0.011 


Table 1. Exact and estimate of the temperature and heat flux 


in the above problem 
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One can see from the data in table 1, the relation errors generated through the 


computation show that the approximate and the exact solutions are vanished. 


Example 4.2. Suppose that 


f(x) = 2° 46a, 
g(t) = 0, 

p(t) = 6t, 

i i | 


Clearly, the exact solution to this problem is u(x,t) = «° + 6at. 


Now, for « = 1, Aty = ay = 0.9, and By = 5, we obtain the following 


a 
30? 
result given in table 2. 


t | at) | u(t) | BGO O! | ua(1,t) | de(,t) | Heo eel 
0 | 0.992943 1 0.007 3 2.9757 0.008 
0.2 | 2.17963 2.2 0.009 4.2 4.13599 0.001 
0.4 | 3.37199 3.4 0.008 5.4 5.31314 0.001 
0.6 | 4.57445 4.6 0.005 6.6 6.52107 0.001 
0.8 | 5.7853 5.8 0.002 7.8 7.75477 0.005 


Table 2. Exact and estimate of the temperature and heat flux 


in Example 4.2 
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